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ABSTRACT 

Several highly luminous Type la supernovae (SNe la) have been discovered. Their high lumi- 
nosities are difficult to explain with the thermonuclear explosions of the Chandrasekhar-mass white 
dwarfs (WDs). In the present study, we estimate the progenitor mass of SN 2009dc, one of the 
extremely luminous SNe la, using the hydrodynamical models as follows. Explosion models of super- 
Chandrasekhar-mass (super-Ch-mass) WDs are constructed, and multi-color light curves (LCs) are 
calculated. The comparison between our calculations and the observations of SN 2009dc suggests 
that the exploding WD has a super-Ch mass of 2.2-2.4 M Q , producing 1.2-1.4 M of 56 Ni, if the 
extinction by its host galaxy is negligible. If the extinction is significant, the exploding WD is as 
massive as ^2.8 M@, and ~1.8 M Q of 56 Ni is necessary to account for the observations. Whether the 
host- galaxy extinction is significant or not, the progenitor WD must have a thick carbon-oxygen layer 
in the outermost zone (20-30% of the WD mass), which explains the observed low expansion velocity 
of the ejecta and the presence of carbon. Our estimate on the mass of the progenitor WD, especially 
for the extinction-corrected case, is challenging to the current scenarios of SNe la. Implications on 
the progenitor scenarios are also discussed. 

Subject headings: supernovae: individual (SN 2009dc) — radiative transfer — white dwarfs 



1. INTRODUCTION 

It has been widely accepted that a Type la super- 
nova (SN la) results from a thermonuclear explosion of 
a carbon-oxygen (C+O) white dwarf (WD) in a close 
binary system. The most likely model is that the explo- 
sion is triggered by carbon ignition in the central region 
of the WD when the WD mass (Mwd) reaches the crit- 
ical mass (Mi a , ~1.38 M© for a non-rotating C+O WD; 
e.g., Hillebrandt & Niemeyer 2000; Nomoto et al. 1997, 
2000). Since Mr a is very close to the Chandrasekhar's 
limiting mass (Chandrasekhar mass, Met) 8 , the result- 
ing explosions are expected to have similar properties. 
Actually, normal SNe la are used as standard candles 
in cosmology (Riess et al. 1998; Perlmutter et al. 1999), 
after correcting their luminosity dispersion by using the 
Pskovskii-Phillips relation (e.g., Pskovskii 1977; Phillips 
1993). 

Despite their uniformity, several unusual SNe la have 
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been found to be much more luminous than normal ones. 
They are SN 2003fg (Howell et al. 2006), SN 2006gz 
(Hicken et al. 2007), SN 2007if (Scalzo et al. 2010; Yuan 
et al. 2010), and SN 2009dc (Yamanaka et al. 2009; 
Tanaka et al. 2010; Silverman et al. 2011; Taubcnberger 
et al. 2011). These SNe la all show slow luminosity 
evolutions (e.g., Scalzo et al. 2010, Table 4). Three 
of them, except for SN 2003fg, show the clear absorp- 
tion line of C II in their early spectra, which are rarely 
detected for normal SNe la (e.g., Marion et al. 2006; 
Tanaka et al. 2008). Such extremely high luminosities 
require >1.2 M Q of radioactive 56 Ni if their explosions 
are spherically symmetric. In order to produce such a 
large amount of 56 Ni, their progenitor C+O WDs are 
suggested to have super- Chandrasekhar (super-Ch) mass 
(i.e., Mwd > Mch), because the exploding WDs should 
contain more than ~0.3 M of the Si-rich layer and the 
unburned C+O layer on top of the 56 Ni-rich core (Howell 
et al. 2006; Hicken et al. 2007; Scalzo et al. 2010; Yuan 
et al. 2010; Yamanaka et al. 2009; Silverman et al. 2011; 
Taubenberger et al. 2011). Alternatively, it could also 
be possible to explain the extremely luminous SNe la 
by asymmetric explosions of Chandrasekhar-mass C+O 
WDs (Hillebrandt et al. 2007). In this paper, we focus on 
SN 2009dc and approximate it with a spherically sym- 
metric model, because the spectropolarimetric observa- 
tions of SN 2009dc suggest that it is a globally spherical 
explosion (Tanaka et al. 2010). 

A super-Ch-mass C+O WD model can be formed if 
it is supported by rapid rotation. For example, Hachisu 
(1986) constructed two-dimensional models of rapidly ro- 
tating WDs. Uenishi et al. (2003) calculated the struc- 
ture and evolution of two-dimensional C+O WDs that 
rotate by getting angular momentum form accreting mat- 
ter. Yoon & Langer (2005) investigated the stability of 
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rapidly rotating C+O WDs for a wider parameter range. 

Explosions and nucleosynthsis of super-Ch-mass C+0 
WDs were simulated by Steinmetz ct al. (1992), Pfannes 
et al. (2010a), and Pfannes et al. (2010b). Maeda 
& Iwamoto (2009) studied the bolometric light curves 
(LCs) of super-Ch-mass WD models. They constructed 
the homologously expanding models of super-Ch-mass 
WDs with parameters of WD mass, 56 Ni mass, abun- 
dance distribution, and so on. Scalzo et al. (2010) studied 
the properties of SN 2007if, assuming a shell-surrounded 
super-Ch-mass WD model as a result of a WD merger. 
They estimated that the ejecta mass is 2.4 M with 1.6 
M Q of 56 Ni. 

By applying Arnett's law to the synthesized bolomet- 
ric LCs, the Ni mass of a SN can be estimated (e.g., 
Arnett 1982). Yamanaka et al. (2009) and Silverman 
et al. (2011) have suggested that SN 2009dc has ~1.2 
M Q if they neglect the extinction by its host galaxy. 
They have also reported that the 56 Ni mass could be as 
large as ~1.6— 1.7 Mq by taking the extinction into ac- 
count. Taubenbcrger et al. (2011) has also analytically 
estimated that the total mass of SN 2009dc is ~2.8 M 
and that the ejected 56 Ni mass is ~1.8 Mq. These ex- 
treme values challenge current models and scenarios for 
SNe la. 

As described above, all the past works rely on the bolo- 
metric LCs, which involves some uncertainties (see Sec- 
tion 2.2). In most cases, the analytic method is used to 
estimate the masses of ejecta and ejected 56 Ni. To derive 
more accurate properties of the super-Ch candidates for 
discussing their progenitor scenarios, more sophisticated 
models are needed. In this paper, we calculate multi- 
color LCs for homologously expanding models of super- 
Ch-mass WDs for the first time. Section 2 describes our 
super-Ch-mass WD models and LC calculations. In Sec- 
tion 3, we compare our results with the photometric and 
spectroscopic observations of SN 2009dc to estimate the 
masses of the ejecta and 56 Ni. Implications of our results 
are discussed in Section 4. We summarize our conclusions 
in Section 5. 

2. MODELS AND CALCULATIONS 

In this section, we describe the procedures for con- 
structing our super-Ch-mass WD models and the code 
for calculating their LCs, to compare with the observa- 
tions of SN 2009dc. Since SN 2009dc has a continuum 
polarization as small as the normal SNe la (Tanaka et al. 
2010), spherical symmetry is assumed. 

2.1. Super- Chandrasekhar- Mass White Dwarf Models 

To construct the homologously expanding models of 
super-Ch-mass WDs, we apply the approximations sim- 
ilar to those adopted by Maeda & Iwamoto (2009). The 
models are described by the following parameters. 

• Mwd- total WD mass. 

• M 5(i Ni (or /56Ni): mass (or mass fraction) of 56 Ni. 
The mass fraction hereafter means the ratio to 

M\VD; e.g., /56Ni = Af56Ni/A^WD- 

• /ece: mass fraction of electron-captured elements 
(ECEs; mostly 54 Fe, 56 Fe, 55 Co, and 58 Ni). Stable 
Fe, Co, and Ni are included, all of which are simply 
assumed to have the same mass fraction. 



/ime ; mass fraction of intermediate- mass elements 
(IMEs). Si, S, and Ca are included, whose mass 
fraction ratio is 0.68 : 0.29 : 0.03, similar to that in 
the Chandrasekhar-mass WD model, W7 (Nomoto 
ct al. 1984; Thielemann et al. 1986). 

fco'- mass fraction of C+O. C and O are contained 
equally in mass fraction. 



Since the equation 



/ECE + /56Ni + /] 



IME 



fee 



(1) 



holds, we eliminate /ime in the following discussion. In 
total, we have four independent parameters. 

The expansion model of a super-Ch-mass WD is con- 
structed as follows. Firstly, with the above parameters, 
we calculate the nuclear energy release during the explo- 
sion (E DUC ) by 



E n 



10 51 erg 



= [0.5/ece + 0.32/ 56N i + 1-24 (1 - / co )] 



(e.g., Maeda & Iwamoto 2009). Then, the binding energy 
of the WD (.Ebin) is evaluated by Equations (22) and 
(32)-(34) described in Yoon & Langer (2005). Here, p c — 
3 x 10 9 g cm~ 3 for all models as in Maeda & Iwamoto 
(2009), so that i?bin increases almost linearly with Mwd- 
For Mwd > 2.1 M Q , we extrapolate the formula of E^- ln 
(Jeffery et al. 2006). These E nuc and -Ebin give the kinetic 
energy of the exploding WD as 



-Ekin — E n 



Ebin- 



(3) 



How much mass fraction of 56 Ni is synthesized at the 
deflagration or detonation wave depends mainly on the 
temperature and thus on the density at the flame front. 
For W7, e.g., 56 Ni is synthesized at the flame densities 
of p <~ p c -2 x 10 7 g cm~ 3 . Then the total synthesized 
56 Ni mass depends on the mass contained in this density 
range, and thus on the presupernova density structure 
of the WD and the flame speed. The rotating WD is 
more massive than the non-rotating WD with the same 
central density (Yoon & Langer 2005). For the same 
central density, therefore, the density profile is shallower 
for more massive stars, and thus the mass contained in 
the density range for 56 Ni synthesis is larger; i.e., a more 
massive WD tends to synthesize more 56 Ni. Also, for the 
same central density, faster flame produces more 56 Ni be- 
cause of less pre-expansion, and the flame speed depends 
on the WD mass, rotation law, density structure, and 
possible transition to detonation. To provide constraints 
on these physical processes, we calculate the light curve 
models with various nuclear yields and -Ekin for the same 
central density models. 

Next, the structure of the explosion model of a super- 
Ch-mass WD is obtained by scaling W7, the canonical 
Chandrasekhar-mass WD one, in a self-similar way. We 
scale the density (p) and velocity (v) of each radial grid 
in the model as 

(4) 



.,5/2 n-3/2 



and 



.,-1/2 j-a/2 

v oc M WD 7 E k < n . 



(5) 
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We then consider the distribution of four element 
groups (ECEs, 56 Ni, IMEs, and C+O). We assume that 
56 Ni in a super-Ch-mass WD model is mixed. The inner 
boundary of the mixing region is set to be v = 5000 km 
s _1 for all the models, because the observed line velocities 
of Si II are >5000 km s _1 (Yamanaka et al. 2009). And 
we set the outer boundary of the mixing region based on 
the mass coordinate; M r = 1.13 (Mwo/Ma) Mq. Here 
the reference mass coordinate (M r = 1.13 Mq) is the 
outer boundary of the 56 Ni distribution in the W7 model. 
This outer boundary corresponds to that of the 56 Ni- 
produced zone in the W7 model. By these definitions, the 
mixing region is uniquely set when we choose four model 
parameters. Note that, for our models, the velocity at 
the outer boundary of the mixing region differs among 
models even with the same Mwd, due to the scaling by 
Equation (5). This mixing is important to account for 
the observed velocities of the Si II line (see Section 3.2). 

To demonstrate the scaling and mixing, we plot the 
structure and abundance distribution of the two super- 
Ch-mass WD models in Figure 1. The model plotted on 
the left panels had no IMEs at 5000 < v < 11000 km 
s _1 before being mixed; the mixing has extended 56 Ni 
outwards and IMEs inwards. In some models, C+O are 
also mixed inwards as in the model right in Figure 1. 

We set the parameter range as follows to compare the 
models with the observations of SN 2009dc. 

• M WD /M Q = 1.8, 2.0, . . ., 2.8. 

• M 56Ni /M Q = 1.2, 1.4, . . ., 1.8, 1.9, 2.0. 

• /ece - 0.1, 0.2, 0.3. 

• /co -0.1, 0.2, 0.3, 0.4. 

Here, /ime is obtained by using Equation (1). Any mod- 
els which have /ime < are not constructed (e.g., models 
with M WD = 2.0 M , /ece = 0.2, and / co > 0.2). 
We also excluded a model where the inner boundary 
of the IMEs layer without mixing is located outer than 
the mixing region (e.g., a model with Mwd = 2.2 M , 
/ece = 0.3, and /co = 0.1), because no IMEs extend in- 
wards by mixing. The parameters and i?kin of our models 
are listed in Columns 2-5 of Table 1 . 

2.2. Light Curve Calculations 

Multi-color LCs for the constructed models are calcu- 
lated with STELLA code (Blinnikov et al. 1998, 2000, 
2006), which solves the one-dimensional equations of ra- 
diation hydrodynamics. STELLA was first developed to 
calculate LCs of Type ILL supernovac (Blinnikov & Bar- 
tunov 1993). Its applications to SNe la are described 
in Blinnikov & Sorokina (2000), Blinnikov et al. (2006), 
and Woosley et al. (2007). In the present study, we 
use homologously expanding ejecta as input (see Sec- 
tion 2). Then the radiation hydrodynamics calculation 
is performed for each model. Maeda & Iwamoto (2009) 
used the one-dimensional gray radiation transfer code 
(Iwamoto 1997; Iwamoto et al. 2000), where they as- 
sumed simplified opacity for line scatterings. STELLA, 
on the other hand, considers 155000 spectral lines in LTE 
assumption to calculate the opacity with expansion effect 
more realistically, as well as free-free/bound-free transi- 
tions and Thomson scattering. 



We note here that the "bolometric" luminosity re- 
ported from the observations is the uvoir luminosity 
{.Luvoir) • Yamanaka et al. (2009) derived Luvoir 

by as- 
suming that (actually observed) integrated BVRL lumi- 
nosity (Lbvri) is 60% of L uvoir . This is commonly ap- 
plied for normal SNe la (Wang et al. 2009a). Since it 
is still unknown whether this assumption is applicable 
to super-Ch candidates, we directly compare theoretical 
and observed Lbvri rather than L UVO i r . This is advan- 
tage of our multi-color calculations. In the calculations, 
the bolometric, uvoir, and BVRL LCs cover 1-50000 A, 
1650-23000 A, and 3850-8900 A, respectively. 

2.3. Extinction 

To compare the calculated and observed LCs, the ob- 
servational data must be corrected for the extinction. In 
order to correct the extinction, three values are needed; 
the B — V excesses for the Milky Way (E(B-V)mw) and 
the host galaxy of SN 2009dc (E(B -V) host ), and the ra- 
tio of the F-band extinction to the B — V excess (Rv) of 
the host galaxy (Rv is set to be 3.1 for the Milky Way). 
Of these values, E(B — V)\ los t and Rv of the host galaxy 
are somewhat difficult to estimate. 

The observed Na I absorption line in the host galaxy 
suggests that the reddening caused by it is not negli- 
gible. To estimate E(B — V)h os t, one may use the ob- 
served color. However, the observed B — V of SN 2009dc 
is significantly different from normal SNe la, which may 
suggest that the Lira-Phillips relation (Lira 1996) should 
not be applied. The other method to derive E(B ~V)h os t 
is using the equivalent width of the Na I D absorption 
line (EWNaiD, e.g. Turatto et al. 2003). But it is also 
noted that the observed EWNaiD has a (relatively) large 
error (Silverman et al. 2011; Taubenberger et al. 2011). 
Also for Rv, a non-standard, smaller value (<3.1) may 
be preferred for normal SNe la (e.g. Nobili & Goodbar 
2008; Wang et al. 2009b; Folatelli et al. 2010; Yasuda & 
Fukugita 2010). If the host galaxy of SN 2009dc also has 
Rv < 3.1, its extinction is overestimated. 

To cover most of possible ranges of extinction, we con- 
sider two extreme cases, where the extinction by the host 
galaxy is negligible (E(B-V)host = mag) and signif- 
icant (E(B-V)host = 0.14 mag), respectively. We set 
E(B — V)mw = 0.71 mag, and Rv = 3.1 for the host 
galaxy (same as for the Milky way). An extinction law 
by Cardelli et al. (1989, Table 3) is applied. 

3. RESULTS 

3.1. BVRL Light Curves 

The left panel in Figure 2 shows the calculated BVRL 
LCs for the super-Ch-mass WD models with different 
Mwd- The other parameters are set to be the same 
(M 56Ni = 1.2 Mq, /ece = 0.1, and / co = 0.2). A clear 
relation is seen between the BVRL LCs and MyvD from 
this panel; a more massive model shows a broader BVRI 
LC. 

Such a mass dependence is consistent with the relation 
between the timescale of the bolometric LC (iboi) and 
M WD , 

t^^R^M^E^J 2 (6) 

(Arnett 1982), where R is the opacity averaged in the 
ejecta (although R does differ among models and with 
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Table 1 

Model Summary 



Name 


M WD 
(M e ) 




/bcb 


fco 


-E-kin 
(10 51 erg) 


WRMSR a 
w/o cxt. b w/ ext. c 


v ph range d 
(10 3 km s- 1 ) 


W7 


1.378° 


0.58 c 


0.22° 


0.14° 


1.3° 


2.51 


3.37 


13.40- 


-7.50 




1.8 


1.2 


0.1 


0.1 


1.49 


0.76 


2.18 


12.58- 


-7.83 




1.8 


1.2 


0.1 


0.2 


1.26 


0.62 


2.10 


12.06- 


-8.34 




2.0 


1.2 


0.1 


0.1 


1.48 


0.51 


2.03 


12.03- 


-8.16 




2.0 


1.2 


0.1 


0.2 


1.23 


0.43 


1.97 


11.37- 


-8.17 




2.0 


1.2 


0.2 


0.1 


1.58 


0.86 


2.23 


12.86- 


-8.14 




2.0 


1.4 


0.1 


0.1 


1.54 


0.41 


1.80 


11.32- 


-7.52 




2.2 


1.2 


0.1 


0.1 


1.46 


0.46 


1.98 


11.53- 


-8.30 




2.2 


1.2 


0.1 


0.2 


1.19 


0.31 


1.92 


10.61- 


-7.82 


A 


2.2 


1.2 


0.1 


0.3 


0.92 


0.42 


1.96 


9.34- 


7.05 




2.2 


1.2 


0.2 


0.1 


1.57 


0.77 


2.17 


12.10- 


-8.45 




2.2 


1.2 


0.2 


0.2 


1.30 


0.62 


2.12 


11.24- 


-8.26 




2.2 


1.4 


0.1 


0.1 


1.53 


0.39 


1.66 


11.56- 


-7.90 




2.2 


1.4 


0.1 


0.2 


1.26 


0.34 


1.61 


10.85- 


-8.09 




2.4 


1.2 


0.1 


0.1 


1.45 


0.32 


1.87 


11.10- 


-8.16 




2.4 


1.2 


0.1 


0.2 


1.15 


0.42 


1.88 


10.05- 


-7.46 


B 


2.4 


1.2 


0.1 


0.3 


0.85 


0.59 


1.99 


8.59- 


6.51 




2.4 


1.2 


0.2 


0.1 


1.57 


0.65 


2.10 


11.68- 


-8.57 




2.4 


1.2 


0.2 


0.2 


1.27 


0.51 


2.07 


10.74- 


-8.05 




2.4 


1.2 


0.3 


0.1 


1.69 


0.95 


2.29 


12.27- 


-8.58 




2.4 


1.4 


0.1 


0.1 


1.51 


0.47 


1.55 


11.22- 


-8.05 




2.4 


1.4 


0.1 


0.2 


1.22 


0.42 


1.56 


10.25- 


-7.72 


C 


2.4 


1.4 


0.1 


0.3 


0.92 


0.52 


1.63 


8.91- 


6.99 




2.4 


1.4 


0.2 


0.1 


1.63 


0.49 


1.79 


11.83- 


-8.15 




2.4 


1.4 


0.2 


0.2 


1.34 


0.31 


1.76 


10.90- 


-8.19 




2.4 


1.6 


0.1 


0.1 


1.58 


0.95 


1.21 


11.50- 


-7.62 




2.4 


1.6 


0.1 


0.2 


1.28 


0.90 


1.24 


10.44- 


-7.90 




2.6 


1.2 


0.1 


0.1 


1.43 


0.35 


1.88 


10.61- 


-7.84 


D 


2.6 


1.2 


0.1 


0.2 


1.11 


0.67 


1.86 


9.54- 


7.01 


E 


2.6 


1.2 


0.1 


0.3 


0.79 


0.86 


2.22 


7.64- 


5.11 




2.6 


1.2 


0.1 


0.4 


0.47 


1.61 


2.74 


5.92- 


4.19 




2.6 


1.2 


0.2 


0.1 


1.56 


0.46 


2.03 


11.20- 


-8.43 




2.6 


1.2 


0.2 


0.2 


1.24 


0.40 


1.99 


10.22- 


-7.69 


F 


2.6 


1.2 


0.2 


0.3 


0.92 


0.52 


2.07 


8.68- 


6.78 




2.6 


1.4 


0.1 


0.1 


1.50 


0.49 


1.52 


10.71- 


-8.06 


G 


2.6 


1.4 


0.1 


0.2 


1.17 


0.56 


1.52 


9.60- 


7.35 


H 


2.6 


1.4 


0.1 


0.3 


0.85 


0.73 


1.64 


8.19- 


6.42 




2.6 


1.4 


0.2 


0.1 


1.63 


0.51 


1.75 


10.98- 


-8.36 




2.6 


1.4 


0.2 


0.2 


1.30 


0.36 


1.65 


10.47- 


-7.99 




2.6 


1.6 


0.1 


0.1 


1.56 


1.05 


1.15 


10.89- 


-7.79 




2.6 


1.6 


0.1 


0.2 


1.24 


1.35 


1.24 


10.21- 


-7.73 




2.6 


1.6 


0.2 


0.1 


1.69 


0.70 


1.41 


11.59 


-7.82 




2.6 


1.8 


0.1 


0.1 


1.62 


1.41 


0.91 


11.01- 


-7.40 




2.6 


1.8 


0.1 


0.2 


1.30 


1.44 


0.94 


10.16- 


-7.69 




2.8 


1.2 


0.1 


0.1 


1.41 


0.55 


1.83 


10.21- 


-7.50 


I 


2.8 


1.2 


0.1 


0.2 


1.07 


0.80 


1.89 


8.98- 


6.30 




2.8 


1.2 


0.1 


0.3 


0.72 


1.48 


2.70 


6.40- 


4.52 




2.8 


1.2 


0.1 


0.4 


0.37 


1.90 


2.94 


4.97- 


3.66 




2.8 


1.2 


0.2 


0.1 


1.55 


0.39 


1.99 


10.88- 


-8.23 


J 


2.8 


1.2 


0.2 


0.2 


1.21 


0.40 


1.97 


9.69- 


7.40 


K 


2.8 


1.2 


0.2 


0.3 


0.86 


0.59 


2.06 


8.09- 


6.26 




2.8 


1.2 


0.3 


0.1 


1.69 


0.73 


2.17 


11.51- 


-8.64 




2.8 


1.2 


0.3 


0.2 


1.35 


0.63 


2.16 


10.35- 


-7.97 




2.8 


1.4 


0.1 


0.1 


1.48 


0.65 


1.48 


10.25- 


-7.85 


L 


2.8 


1.4 


0.1 


0.2 


1.13 


0.81 


1.52 


9.18- 


6.98 


M 


2.8 


1.4 


0.1 


0.3 


0.78 


u.yi 


liIO 


T KK 
1 .OD- 


o.oo 




2.8 


1.4 


0.2 


0.1 


1.62 


0.33 


1.66 


10.98 


-8.27 




2.8 


1.4 


0.2 


0.2 


1.27 


0.58 


1.62 


10.08- 


-7.79 




2.8 


1.4 


0.3 


0.1 


1.76 


0.52 


1.86 


11.57- 


-8.46 




2.8 


1.6 


0.1 


0.1 


1.54 


1.13 


1.13 


10.42- 


-7.81 




2.8 


1.6 


0.1 


0.2 


1.20 


1.34 


1.18 


9.54- 


7.37 




2.8 


1.6 


0.1 


0.3 


0.85 


1.17 


1.37 


7.90- 


6.35 




2.8 


1.6 


0.2 


0.1 


1.68 


0.74 


1.35 


11.15- 


-8.00 




2.8 


1.6 


0.2 


0.2 


1.34 


0.80 


1.31 


10.14- 


-7.92 




2.8 


1.8 


0.1 


0.1 


1.61 


1.68 


0.80 


10.77- 


-7.61 


N 


2.8 


1.8 


0.1 


0.2 


1.26 


1.70 


0.88 


9.72- 


7.57 




2.8 


1.9 


0.1 


0.1 


1.64 


1.78 


0.74 


10.68- 


-7.45 





2.8 


1.9 


0.1 


0.2 


1.29 


1.85 


0.76 


9.80- 


7.60 




2.8 


2.0 


0.1 


0.1 


1.67 


1.98 


0.62 


10.64- 


-7.24 



Note. — The parameters and values of the selected models are bold-faced. 
a Weighted root-mean-square residual of the BVRI LC, calculated by Equation (7). 
b The extinction by the host galaxy of SN 2009dc is neglected (E(B— V)host = mag). 
c The host-galaxy extinction is corrected (E(B — V)host = 0.14 mag). 
d Maximum and minimum values during —5 < t_e < 25 days. 
c Nomoto et al. (1984); Thielemann et al. (1986). 
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Figure 1. Model configurations of two super-Ch-mass WD models as examples. The parameters of these models are A/wd = 2.0 Mq, 
M 56Ni = 1.2 M Q , /ece = 0.2, and / C o = 0.1 in the left panels, and M W D = 2.4 M Q , M 56N i = 1.2 M Q , /ece = 0.1, and / C o = 0.3 
in the right panels. Top: density distribution of the model (red) and W7 (Chandrasekhar-mass WD model; black) at 2 X 10 4 s after the 
explosion. Bottom: abundance distribution of the model. The lines indicate ECEs (blue), 56 Ni (green), IMEs (orange), and C+O (black). 
The mixing region is expressed as the shaded area. 



I 





20 40 60 

days since explosion 

Figure 2. Left: BVRI LCs of the super-Ch-mass WD models are plotted (solid lines). The parameters of the plotted models are 
M 56 Ni = 1-2 M Q , / ece = 0.1, and f co = 0.2. M WD of the models are 1.8 M (red), 2.0 M (blue), 2.2 Mq (green), 2.4 Mq (orange), 
2.6 Mq (magenta), and 2.8 Mq (grey). The BVRI LC of W7 is also plotted (black dotted). Right: lines are same as the left plot, but the 
x-axis is now ts- In addition, the observed LC of SN 2009dc is also shown. The squares are the observations by Yamanaka et al. (2009), 
where the error bars are set to be 20% of LgvRl- The triangle around tg ~ —20 days indicates the early detection in Silverman et al. 
(2011), assuming that the ij-band luminosity is 20% of the L uvo i r (equivalent to LgvRl /3) and E(B — V)mw = 0.07 mag. 
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time). Note that E^ n anticorrelates with Mwd for the 
models plotted in the left panel of Figure 2 because /ece 
and fco are fixed (c.f. Table 1). 

We consider a timescale to see quantitatively if the 
BVRI LCs depends on M WD and E^ in , analogous to 
Equation (6). For this purpose, we take the declining 



timescale (t 



+ 1/2) 



which is defined as the time for the 



BVRI LC to halve its luminosity after the peak. In the 
left panels of Figure 3, shown are the MwD~-E : km (large 
panel), MwD^i+i/2 (small, top), and i?km-i+i/2 (small, 
right) plots for the models listed in Table 1. The BVRI 
LCs of the massive and less energetic models tend to be 
wider (i.e. larger t +1 / 2 ), which is expected from the left 
panel in Figure 2. For the dependence of i+1/2 on Mwd 

and -Ekim we plot t+1/2 against M^^E kin ^ 2 in the right 
panel of Figure 3. An almost linear relation is seen be- 
tween t + i/2 and M^Jj^E^ 2 . We thus confirm a relation 
similar to Equation (6) for the BVRI LCs. 

The calculated and observed BVRI LC are shown on 
the right panel in Figure 2, by setting the dates of max- 
imum in the B-band of the models and observations at 
the same day as ts = day. We plot the observational 
data provided by Yamanaka et al. (2009), neglecting the 
extinction by the host galaxy (open squares). The open 
triangle corresponds to the early i?-band detection by Sil- 
verman et al. (2011), where we simply assume that the 
i?-band luminosity is equivalent to 20% of L uvoir , i.e., 
1/3 of Lbvri, assuming L uvoir = L BVRI /0.6 (cf. Wang 
et al. 2009a, Figure 24). The BVRI LCs of the super-Ch- 
mass WD models on the right panel are as luminous as 
the observations around the maximum, while some have 
relatively broader BVRI LCs or shorter rising time than 
the observations. Especially, the two less massive models 
(Mwd = 1.8-2.0 M©) are not preferred because they are 
too faint at ts ~ —20 days. 

In order to find the well-fitted models, we calculate 
the weighted root-mean-square residual (WRMSR) of the 
BVRI LC, 



WRMSR 




(calc) 
BVRI A 



- (obs) 
J BVRI,, 



r r (obs) 
ulj BVRI, 



(7) 



for each model. We use the observational data for 
— 10 < ts < 40 days taken from Yamanaka et al. (2009), 
where the total number of the observations, N, is 23 for 
this time range. The early detection by Silverman et 
al. (2011) is excluded in calculating the WRMSR. Here, 

^bvri i ano - ^bvri i respectively denote L BVRI of the 
observation and calculation obtained at the i-th epoch. 
^^bvri i 1S observational error of L^y^j i , which is 
set to 20% of L { ^l I t . Maeda & Iwamoto (2009) used 
the decline rate of the bolometric LC after maximum 
for comparison between their models and observations. 
We use the above WRMSR instead of the decline rate, 
so that brightness can also be taken into account. The 
WRMSR is calculated with the observational error, and 
shows which model fits relatively to the observed BVRI 
LC. 

We list the WRMSR of the models for the two ex- 
tinction cases in Columns 7 and 8 of Table 1. Good 



agreement is found between the calculated and observed 
BVRI LCs for whole range of Mwd in this study, if the 
extinction by the host galaxy is not considered. In the 
case where the extinction is corrected, most models have 
larger WRMSR (due to smaller Lbvri), while some mod- 
els with large Mwd and M 56 Ni are better fitted. 

3.2. Photospheric Velocity 

Next, we compare the photospheric velocity (v p h) of 
the models and the observed line velocity of Si II. In 
the last column of Table 1, listed are the t> p h ranges of 
our models at — 5 < ts < 25 days. This period covers 
the whole phases of the spectroscopic observations by 
Yamanaka et al. (2009). We estimate the position of the 
photosphere (hence, v p h) in the model from the optical 
depth at the i?-band, where the rest-frame wavelength 
of the observed Si II line (6355 A) is located. 

The velocity of the observed Si II line is reported as 
<9000 km s _1 at the period (Yamanaka et al. 2009). 
Since the Si II absorption line is formed by the Si above 
the photosphere, the calculated v p h should be smaller 
than the observed Si II line velocity (for further details, 
see Tanaka et al. 2011, Figure 1). For the models with 
Mwd = 1-8 Mq and 2.0 Mq, however, the calculated 
v p h are > 10000 km s _1 , much larger than the observed 
line velocity of Si II. Thus, these less massive models are 
far from consistent with SN 2009dc. On the other hand, 
many models with Mwd > 2.2 M Q have v p h < 10000 
km s _1 during — 5<is<25 days, which are somewhat 
compatible with observations. These models with small 
Vph commonly have such a large C+O mass fraction as 
/ co = 0.2-0.4. From Equation (2), E nuc and thus E^ n 
are affected mainly by fco rather than /ece- A model 
with larger fco has smaller E nuc , thus smaller £^in and 
Vph, being preferred for SN 2009dc. 

3.3. Plausible Models for SN 2009dc 

In order to find the most plausible model for SN 
2009dc, we first select models based on BVRI LC and 
Vph discussed above. We use criteria of WRMSR < 1 
and Vph < 10000 km s _1 during the —5 < ts < 25 days. 
By these criteria, models with a label A-0 in Column 1 
of Table 1 are selected. Among them, models A-M are 
selected for the case where the host-galaxy extinction is 
neglected; only models N and for the significant extinc- 
tion. In Figure 4, we plot WRMSR and maximum v p h 
during the period, where the dotted lines indicates the 
criteria. 

We further analyze multi-color LCs of these mod- 
els. Figures 5-9 show the BVRI (left), v p h (right top) 
monochromatic LCs (right bottom) of these models. In 
these figures, the BVRI, B-, V-, R-, and /-band LCs by 
Yamanaka et al. (2009) and the early detection and U- 
band LC by Silverman et al. (2011) are plotted as squares 
and triangles, respectively. Filled and open symbols show 
the observed LCs with and without the extinction cor- 
rection. 

With Figures 5-9, we make more detailed comparisons 
of models A-0 with the observational data of SN 2009dc. 
The comparisons are summarized as follows: 

• The early detection in Silverman et al. (2011) en- 
ables us to constrain the rising time of the models. 
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Figure 4. Plots for the WRMSR calculated by Equation (7) and the maximum v p ^ during —5 < te < 25 days of the models. The crosses 
with the same color correspond to the models with the same Mwd- The two dotted lines indicate the criteria for the model selection (see 
the text). Left: E(B-V) host — mag is assumed; i.e. the extinction by the host galaxy is not considered in calculating the WRMSR. 
Right: the significant extinction (E(B — V)host = 0.14 mag) is assumed. 



The i?-band LCs of models E, H, and M are too 
bright at ts ~ —20 days regardless of whether the 
host-galaxy extinction is corrected or not. These 
models can be excluded because their LCs evolve 
too slowly. 

• The remaining models for the case of the neglected 
host-galaxy extinction (A-D, F, G, and I-L) have as 
small MsgNi as 1.2 M & and 1.4 Mq, which is con- 
sistent with the analytical estimates. While model 
J has the smallest WRMSR, its v p h is relatively 
larger than the observations. Models D, G, and L 
also have larger v p h around ts ~ day. Among the 



other models (A-C, F, I, and K), the models with 
Mwd < 2.4 M Q (A-C) are preferred, because their 
BVRI LCs reproduce the observations well around 
the peak. Especially, model B has the smallest v p h 
of the three models, so we suggest that model B is 
the most plausible model for SN 2009dc without 
the host-galaxy extinction. 

• If the extinction by the host galaxy is significant, 
we have only two candidate models with JWwd = 
2.8 Mq (N and 0). They both have similar proper- 
ties; slightly larger Lbvri at is ~ —20 days, but 
a bit smaller Lbvri around ts ~ day. We re- 
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t B [day] t„ [day] 



Figure 5. BVRI and monochromatic LCs, and i> p h for models A, B, and C. Left: same as the right panel of Figure 2. For the open 
points, we consider the extinction by the host galaxy; for the solid ones, we neglect it. The point at the left-top corner indicates the error 
bars for squares. Right top: with the line velocity of C II (pluses) and Si II (crosses) observed by Yamanaka et al. (2009). Right 
bottom: same as the left panel, but multi-band LCs are plotted. The observation data in (7-band (cyan triangles) and the early detection 
in R-band (red triangles around tg ~ —20 days) are taken from Silverman et al. (2011). Note that we assume E(B — V)mw = 0.07 mag 
and E(B — V)host = (open triangles), 0.10 mag (filled triangles). 
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Figure 6. Same as Figure 5, but for models D, E, and F. 



gard model N as the most plausible model for the 
extinction-corrected case because the rising time is 
shorter than model 0. 

By these results, we suggest that the mass of progen- 
itor WD of SN 2009dc is 2.2-2.4 M if the extinction 
by its host galaxy is negligible, and ~2.8 Mq with the 
extinction if the extinction is significant. The 56 Ni mass 
needed for the SN is 1.2-1.4 Mq for the former case, and 
1.8-1.9 Mq for the latter, respectively. For the latter 
case with extinction, the estimated masses are consis- 
tent with those suggested by Taubenberger et al. (2011), 
though they assume that the mean optical opacity of SN 
2009dc is similar to (normal) SN 2003du. 



For SN 2007if, one of the super-Ch candidates, Scalzo 
et al. (2010) estimate the total mass of its progenitor to 
be ~2.4 Mq, as massive as we estimate for SN 2009dc. 
They consider a shell-structured model to explain the low 
Si II line velocity and its plateau-like evolution, where 
the massive envelope decelerates the outer layers of the 
ejecta. Our calculations also reproduce the lower Si II 
line velocity and similar evolution of SN 2009dc (the 
right-top panels of Figures 5-9), and suggest that the 
low line velocity and flat evolution of Si II can be ex- 
plained by scaled super-Ch-mass WD models, as well as 
the shell-shrouding models. 
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4. DISCUSSION 

4.1. Light Curves of Super-Chandrasekhar-Mass White 
Dwarf Models 

4.1.1. Multi-Color Light Curves 

We have calculated the multi-color LCs of the super- 
Ch-mass WD model for the first time (right-bottom pan- 
els in Figures 5-9). Even the model being in good agree- 
ment with Lbvri and w p h shows discrepancy for each 
band to some extent, especially in the / band. This de- 
viation from observed SNe la in the / band is also seen 
for a Chandrasekhar-mass WD model calculated by the 
code STELLA (Woosley et al. 2007, Figure 9). It could 
be improved by taking more spectral lines into account 
for the opacity calculation. 

The comparison with the velocity suggests that the 
mixing occurred in the ejecta. It is interesting to note 



that LCs of SN 2009dc in the L band do not clearly 
show double-peak features. For Chandrasekhar-mass 
WD models, Kasen (2006) calculated the multi-color 
LCs, to conclude that mixing in the ejecta could produce 
the single-peak LCs in the I band. In fact, our models 
are partially mixed and their /-band LCs do not show 
two peaks (Figures 5-9). It could be also the case for 
super-Ch-mass WD models that mixing affects their I- 
band LCs, although some uncertainties mentioned above 
should be taken into account. 

4.1.2. Bolometric, uvoir, and BVRI Light Curves 

In Figure 10, the bolometric, uvoir, and BVRI LCs are 
plotted for models W7, B, and N. The uvoir LCs peaks 
earlier than the BVRI LCs, by ^2 days for W7 model, 
and by ~10 days for model B and N, respectively. As 
for the peak luminosity, the uvoir LC of W7 is brighter 
than the BVRI LCs by ^0.2 dex, while those of the 
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Figure 9. Same as Figure 5, but for models M, N, and 0. 

two super-Ch-mass WD models, by ~0.3 dex. These 
shifts are understood by considering the ultraviolet (UV) 
radiation in the early phase (Blinnikov & Sorokina 2000). 

Most of BVRI LCs for the super-Ch-mass WD mod- 
els in this paper seems to be fainter for ts < day and 
slightly brighter after that, than the observations (Fig- 
ures 5-10). This discrepancy could be solved by changing 
the 56 Ni distribution in the model, which powers its ra- 
diation. The more 56 Ni outside could shorten the rising 
time, or make more luminous before the peak, while the 
less 56 Ni inside could make fainter after that. Further 
detailed modeling is needed since we simply assume the 
56 Ni distribution analogous to W7. 

In Figure 10, we add the observation data for ts > 90 
days plotted in Figure 7 of Taubenberger et al. (2011) as 
filled pentagons, which are also corrected for the host- 
galaxy extinction by using their values. The calculated 
bolometric LC fits to the observed LC tail quite well. 
The observed tail lie just on model N, which means that 
the estimate on M56M for the model is consistent. 

4.2. Instabilities of Rotating White Dwarfs 

We calculate -Ebin by just extrapolating the formula, 
valid only for M WD > 2.1 M Yoon & Langer (2005). 
Using their Equations (22) and (29)-(31), we can calcu- 
late the ratio of the rotational energy (T) to (the magni- 
tude of) gravitational one of (W) a super-Ch-mass WD 
with our parameter p c . This ratio T/W indicates the 
stability of the WD. If T/W is small enough, the WD 
is stable against rotation. If T/W reaches 0.14, the WD 
suffers the non-axisymmetric instability (e.g. Ostriker & 
Bodenheimer 1973). The above equations show that a 
WD with M W d ~ 2.4 M Q has T/W = 0.14 for our p c . 
A super-Ch-mass WD with Mwd > 2.4 M Q might not 
form. However, T/W as a function of Mwd differs with 
the rotation law (e.g., Hachisu 1986), so that it is pos- 
sible that even such a massive WD as Mwd > 2.4 M 
has T/W < 0.14 and thus form. We thus consider the 
models with A/wd > 2.4Af Q . 

Our results in this paper suggest that the progenitor 
WD mass of SN 2009dc exceeds 2.2 M even for the 
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Figure 10. Calculated bolometric (solid), uvoir (dotted), and 
BVRI (dashed) LCs for models W7 (black), B (red), and N (blue). 
The observed LCs are also shown as the open and filled symbols for 
the cases with and without the host-galaxy extinction, respectively. 
For the open triangle and squares, the squares are same as in the 
right panel of Figure 2. We assume E(B — V)host = 0.10 mag 
for the filled triangles. The pentagons are taken from Figure 7 of 
Taubenberger et al. (2011), where E(B — V)mw = 0.071 mag and 
E(B-V) host =0.10 mag. 



case where the host-galaxy extinction is negligible. Such 
a large Mwd might suggest that the explosion of this 
super-Ch-mass WD is triggered when its mass reaches 
the critical mass for the above instability of the rotating 
WDs (Hachisu et al. 2012). 

4.3. Possible Progenitors and Scenarios 

The progenitor WD mass for SN 2009dc in our estimate 
largely exceeds Mch, thus putting severe constraints on 
the presupernova evolution of the binary system. This 
also should have important implications on the progeni- 
tor scenarios of ordinary SNe la, both the single degen- 
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erate (SD) and double degenerate (DD) scenarios. 

4.3.1. Single Degenerate Scenario 

The initial masses of the WD (Mwd.o) and its com- 
panion star (M comp! o) should be sufficiently large at 
the beginning of accretion in order to increase Mwd 
to >2.4 Mq. Chen & Li (2009) argued that even for 
Mwd,o = 1-24 M© the accretion from the companion 
of M comP;0 < 3.5 M© docs not obtain Mwd ^ 1-7 
Mq. However, they does not take into account the ef- 
fect of mass-stripping from the companion star due to 
the strong WD wind (Hachisu et al. 2007). Because the 
mass-stripping effectively reduces the mass transfer rate 
from the companion to the WD, the companion star can 
be as massive as M com po ~ 4—7 Mq and Mwd can reach 
—2.4 Mq (Hachisu et al. 2012). Still, to realize such a 
massive WD, Mwd,o ^ 1-1 Mq is preferable (Chen & Li 
2009; Hachisu et al.' 2012). 

The formation of the C+O WD with M W d,o > 1-1 
Mq is realized in the special binaries. In stars of main- 
sequence mass of <8 Mq, the C+O core mass is <1.07 
Mq to avoid off-center C-ignition before the AGB phase 
(e.g., Umeda et al. 1999). After the dredge- up of the 
He layer, the C+O core increases its mass during the 
AGB phase, if the binary separation is wide enough to 
accommodate the AGB star. In such a binary system, 
the C+O WD with M W d,o > 1.1 M© can be formed, if 
the C+O core of the AGB star has already grown massive 
when the AGB envelope is lost in a wind or by Roche- 
lobe overflow. Thus Mwd,o is larger if the mass loss 
rate from the AGB star is smaller. Therefore, the C+O 
WDs with larger Mwd,o is more likely to form in the 
lower metallicity system and in the initially wider binary 
(Hachisu et al. 2007). Since the binary must also be close 
enough for the mass-transfer to occur, the suitable binary 
system could be rare, which is consistent with the low 
occurrence frequency of the super-Ch-mass explosion. 

These requirements of large enough Mwd,o and 
M comPj0 in the SD scenario predicts that SNe la from 
super-Ch-mass WDs are associated with the star-forming 
region and low metallicity environment. It is interesting 
to note that the hosts of the observed super-Ch candi- 
dates are faint and star-forming galaxies except for SN 
2009dc (e.g. Taubenberger et al. 2011, Table 7). For SN 
2009dc, the host galaxy UGC 10064 is a passive (SO) 
galaxy. However, at —40 kpc away from UGC 10064, 
there is a blue irregular galaxy UGC 10063, which could 
also have a star formation in the recent past by the in- 
teraction (Silverman et al. 2011). 

4.3.2. Double Degenerate Scenario 

For the DD scenario, to form a WD of >2.4 M© by 
merging of two C+O WDs, the primary C+O WD needs 
to be initially as massive as Mwd,o ^ 1-33 M© because 
of the following reason. To form a massive C+O core in 
the primary AGB star, the initial binary system needs 
to be wide enough. The Roche lobe overflow of the pri- 
mary AGB star is so rapid that a formation of a common 
envelope is unavoidable. After the loss of mass and an- 
gular momentum from the common envelope, a primary 
C+O WD and the secondary star are left in a binary 
with a small separation. If the separation is too small 
for the secondary star to become an AGB star, the mass 
of the secondary C+O WD is M W d,o < 1-07 Mq. In 



order to form a WD of >2.4 Mq, the primary C+O WD 
should be more massive than 1.33 Mq, which would be 
very rare. In the above scenario, formation of the dou- 
ble C+O WDs whose initial masses are both —1.2 M© 
may not be possible because of the shrink of the binary 
system after the first common envelope phase. 

5. CONCLUSIONS 

To constrain the properties of SN 2009dc, we have cal- 
culated multi-band LCs for the exploding super-Ch-mass 
WD models with a range of model parameters. We find 
that the mass of the WD and other model parameters 
are constrained as follows. 

• The observed BVRI LCs of SN 2009dc are well- 
explained by the super-Ch-mass WD models with 
Mwd = 1.8-2.8 M© and M 56N i = 1.2-1.8 M Q , if 
the extinction by the host galaxy is negligible. 

• The observed line velocity of Si II is consistent with 
w p h of several models with Mwd = 2.2-2.8 Mq, 
but significantly lower than v p h of the less massive 
models. 

• Among our models, the most plausible model is 
model B with M W d = 2.4 M© (1.2 M© of 56 Ni, 
0.24 Mq of ECEs, 0.24 M© of IMEs, and 0.72 M© 
of C+O) for the case without the host-galaxy ex- 
tinction. 

• If the extinction is considered, the mass of the 
super-Ch-mass WD needs to be as massive as 
M W D ~ 2.8 M© (i.e. model N with 1.8 M© of 56 Ni, 
0.28 Mq of ECEs, 0.16 M© of IMEs, and 0.56 M© 
of C+O). We find that the fit to the observation is 
less successful for model N than model B. 

• Such a large Mwd might suggest that the explosion 
of the super-Ch- WD might be related to the onset 
of the instability of the differentially rotating WD. 

Our results in this paper are based on the simplified 
models of the super-Ch-mass WDs. There are still sev- 
eral uncertainties in the models and LCs; such as the 
parameterization of the models, the opacity calculated 
by the code, aspherical effects, and effects of possible cir- 
cumstellar interaction, as well as the instability of the 
massive WDs. However, Mwd and M 56 Ni of the plausi- 
ble models for SN 2009dc are quite consistent with the 
observations, suggesting that our present approach works 
well for this super-Ch candidate. 
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